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INTRODUCTION 


This  program  seeks  to  obtain  a  fundamental  understanding  of  the  chemical 
mechanisms  by  which  enzymes  repair  damaged  DNA,  and  to  use  this  information  to 
design  small  molecule  inhibitors  of  these  enzymes.  The  driving  force  for  these  efforts  is 
the  recognition  that  the  effectiveness  of  cancer  chemotherapy  regimes  is  intimately 
connected  to,  and  in  some  cases  directly  relies  on,  DNA  damage  repair  pathways.  A  more 
sophisticated  understanding  of  the  roles  of  DNA  damage  repair  in  the  pharmacology  of 
DNA  replication  inhibitors  will  allow  for  the  design  of  better  treatments  against  breast 
and  other  cancers. 

Studies  during  this  period  have  focused  on  Task  1  and  2  in  the  original  Statement 
of  Work.  We  have  explored  how  uracil  DNA  glycosylase  (UDG)  binds  to  substrate  DNA 
with  such  high  affinity,  and  have  utilized  this  knowledge  to  synthesize  a  very  tight 
binding  DNA-based  inhibitor  of  the  enzyme.  We  also  have  used  a  new  high  throughput 
screen  for  UDG  to  identify  inhibitors  from  a  novel  directed  uracil  library  that  we  have 
constructed.  This  reporting  period  has  produced  one  publication  in  a  peer-reviewed 
journal  and  one  more  manuscript  is  in  preparation. 

BODY 

A  major  accomplishment  of  this  funding  period  was  the  identification  of  the  damaged 
DNA  features  that  are  important  for  high  affinity  binding  by  UDG.  This  work,  which 
was  published  in  the  journal  Biochemistry  (see  Appendix),  revealed  that  tight  binding 
DNA  molecules  could  be  constructed  if  the  uracil  substrate  was  presented  to  the  enzyme 
in  the  context  of  a  base  pair  that  lacks  hydrogen  bonds  and  exhibits  greater  flexibility 
than  normal  DNA  base  pairs.  Such  designer  uracil  sites  bind  43,000-fold  more  tightly 
than  sites  that  lack  these  features. 

Our  new  small  molecule 
inhibitor  development  strategy  uses 
the  specificity  and  binding  energy  of 
the  uracil  ligand  to  target  the  UDG 
active  site,  and  then  relies  upon 
random  exploration  of  binding  sites 
for  tethered  R  groups  (Fig.  1). 

Various  types  of  tethering  strategies 
are  being  pursued  in  the 
pharmaceutical  industry  and 
successful  results  in  several  systems 
have  been  reported.  Our  approach  has 
the  broader  impact  of  exploring  novel 
synthetic  and  screening  strategies  to  discover  inhibitors  of  enzymes  that  repair  DNA. 
Such  targets  are  largely  unexplored  in  terms  of  their  therapeutic  potential. 

This  library  is  constructed  by  attaching  one  end  of  a  bis-(aminooxy)alkane  linker 
scaffold  to  C6,  C5  and  N1  aldehyde  derivatives  of  uracil  as  shown  in  Fig.  2  .  The  uracil- 
oximes  may  then  be  connected  to  a  variety  of  aldehydes  using  the  free  ONH2  group  on 
the  linker.  The  bis-(aminooxy)alkane  scaffold  is  ideal  because  it  allows  variation  of  the 
linker  length  to  probe  different  potential  binding  surfaces  on  the  enzyme  (Fig.  2),  and  the 


UDG 

uracil  site 

uracil  site 

A A 

As 

A 

VA  vfU 

Site  A  Site  B 

Site  A  Site  B 

Figure  1.  Targeting  UDG  using  a  simple  substrate 

tethering  approach. 

Bifunctional  (URR)  and 

monofunctional  (UAR)  ligands  are  depicted. 
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reactions  go  to  completion 
negating  the  requirement 
for  purification  steps  (see 
below).  The  potential  size 
of  the  directed  MF  library 
may  be  estimated  as  ~ 

4,500  compounds.  This 
calculation  assumes  the 
reaction  of  -300  random 
RCHO  compounds  X  3 
uracil  aldehyde  derivatives 
X  5  linker  lengths  =  4,500 
unique  compounds. 

KEY  RESEARCH 

ACCOMPLISHMENTS 

•  Elucidated  the 

molecular 

determinants  for  high  affinity  binding  of  uracil-substrates  by  UDG  and  designed  a 
high  affinity  DNA-based  inhibitor  (SOW,  Task  1). 

•  Designed  and  screened  a  novel  uracil-based  small  molecule  library  for  inhibition 
of  hUDG  (SOW,  Task  2). 

•  Identified  and  characterized  several  small  molecule  inhibitors  of  hUDG. 
REPORTABLE  OUTCOMES: 

This  work  has  resulted  in  one  published  manuscript: 

1 .  Krosky,  D.  J.  and  Stivers,  J.  T.  (2005)  The  Origins  of  High  Affinity  Enzyme  Binding 
to  an  Extrahelical  DNA  Base,  Biochemistry  44,  5949-5959. 

CONCLUSIONS: 

The  long-term  goal  of  this  research  is  to  increase  the  effectiveness  of  5-FU  chemotherapy 
through  the  selective  targeting  of  BER  repair  enzymes  such  as  UDG.  We  have  made 
significant  progress  towards  this  goal  in  this  funding  period  through  the  synthesis  of  a 
very  potent  DNA-based  inhibitor  of  the  enzyme.  We  have  also  made  significant  progress 
towards  the  development  of  cell  permeable  small  molecule  inhibitors  of  UDG. 

REFERENCES:  None 

APPENDICES:  One  manuscript  is  included  (see  above). 
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Figure  2.  HTP  synthesis  of  directed  uracil  library  for  hUDG 
inhibition  using  oxime  chemistry.  (A)  Mixed  oxime  formation 
involves  the  reaction  of  an  aikyldioxyamine  of  variable  length  with  a 
uracil  aldehyde  and  a  RCHO  library.  The  reaction  proceeds 
quantitatively  to  produce  a  statistical  mixture  of  oxime  derivatives  (50 
%  UAR)  .  The  O-alkyl  oxime  libraries  are  rapidly  and  efficiently 
synthesized  in  96  well  plates  without  purification. 


5 


Biochemistry 


©  Copyright  2005  by  (he  American  Chemical  Society 


Volume  44,  Number  16 


April  26,  2005 


Articles 
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abstract:  Base  flipping  is  a  highly  conserved  strategy  used  by  enzymes  to  gain  catalytic  access  to  DNA 
bases  that  would  otherwise  be  sequestered  in  the  duplex  structure.  A  classic  example  is  the  DNA  repair 
enzyme  uracil  DNA  glycosylase  (UDG)  which  recognizes  and  excises  unwanted  uracil  bases  from  DNA 
using  a  flipping  mechanism.  Previous  work  has  suggested  that  enzymatic  base  flipping  begins  with  dynamic 
breathing  motions  of  the  enzyme-bound  DNA  substrate,  and  then,  only  very  late  during  the  reaction 
trajectory  do  strong  specific  interactions  with  the  extrahelical  uracil  occur.  Here  we  report  that  UDG 
kinetically  and  thermodynamically  prefers  substrate  sites  where  the  uracil  is  paired  with  an  unnatural 
adenine  analogue  that  lacks  any  Watson— Crick  hydrogen-bonding  groups.  The  magnitude  of  the  preference 
is  a  striking  43000-fold  as  compared  to  an  adenine  analogue  that  forms  three  H-bonds.  Transient  kinetic 
and  fluorescence  measurements  suggest  that  preferential  recognition  of  uracil  in  the  context  of  a  series  of 
incrementally  destabilized  base  pairs  arises  from  two  distinct  effects:  weak  or  absent  hydrogen  bonding, 
which  thermodynamically  assists  extrusion,  and,  most  importantly,  increased  flexibility  of  the  site  which 
facilitates  DNA  bending  during  base  flipping.  A  coupled,  stepwise  reaction  coordinate  is  implicated  in 
which  DNA  bending  precedes  base  pair  rupture  and  flipping. 


Enzymes  that  modify  or  cleave  nucleobases  in  DNA,  such 
as  DNA  methyltransferases  and  DNA  glycosylases,  are 
confronted  with  a  formidable  chemical  problem:  gaining 
access  to  substrate  bases  that  are  sequestered  inside  the  DNA 
double  helix  (1—3).  A  conserved  enzymatic  solution  to  this 
problem  is  base  flipping,  where  the  target  base  and  sugar 


t  This  work  was  supported  by  National  Institutes  of  Health  Grant 
GM56834-10.  D.J.K.  was  supported  by  the  DOD  Breast  Cancer 
Research  Program  (DAMD1 7-03- 1-1251). 

*  To  whom  correspondence  should  be  addressed.  Tel:  410-502-2758. 
Fax:  410-955-3023.  E-mail:  jstivers@jhmi.edu. 

*  The  Johns  Hopkins  School  of  Medicine. 

5  Center  for  Advanced  Research  in  Biotechnology  of  the  University 
of  Maryland  Biotechnology  Institutes  and  the  National  Institute  of 
Standards  and  Technology. 


are  extruded  from  the  DNA  duplex  into  the  enzyme  active 
site  (4). 

Structural  and  mechanistic  studies  indicate  that  base 
flipping  is  a  multistep  process  involving  two  coupled  reaction 
coordinates  (Figure  1  A)  (5—12),  which  are  depicted  using  a 
free  energy  contour  plot  in  Figure  IB.  The  first  coordinate 
involves  ~180°  rotation  of  the  entire  target  nucleotide  from 
the  duplex  stack  (vertical  axis,  Figure  IB),  while  the  second 
coordinate  involves  enzyme-induced  DNA  bending  (hori¬ 
zontal  axis,  Figure  IB).  The  progress  along  each  reaction 
coordinate  may  be  perfectly  synchronized  (diagonal  dashed 
line,  Figure  IB),  or  alternatively,  one  process  may  proceed 
ahead  of  the  other  (curved  trajectory,  Figure  IB).  For 
instance,  if  DNA  bending  proceeds  ahead  of  nucleotide 
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Figure  1 :  Uracil  flipping  has  two  coupled  reaction  coordinates  involving  base  rotation  and  DNA  bending.  (A)  Upon  binding  to  free  B 
DNA  (left)  UDG  rotates  the  deoxyuridine  nucleotide  from  the  DNA  base  stack  by  180°  and  bends  the  DNA  by  about  40°.  The  DNA 
structure  on  the  right  was  extracted  from  the  complex  of  UDG  with  substrate  analogue  DNA  (PDB  code  1EMH)  (8).  (B)  A  two-dimensional 
free  energy  contour  map  depicting  the  two  coupled  reaction  coordinates  of  base  rotation  (vertical  coordinate)  and  DNA  bending  (horizontal 
coordinate).  The  diagram  shows  the  two  enzyme-bound  intermediates  (1|,  I2)  that  have  been  previously  detected  using  rapid  kinetic  methods 
(5,  6)  and  the  final  product  of  the  flipping  and  bending  reaction  (P).  In  principle,  progress  along  both  reaction  coordinates  may  be  synchronized 
(diagonal  line),  or  one  process  may  lag  behind  the  other.  The  contour  map  depicts  the  reaction  trajectory  (curved  dashed  line)  where 
bending  precedes  flipping,  and  a  low  energy  (blue)  bent  intermediate  is  formed  before  significant  progress  along  the  base  rotation  coordinate 
occurs  (I2).  Thus,  flipping  uracil  from  unbent  DNA  is  a  high-energy  (red)  improbable  process  (upper  left  comer).  Formation  of  I2  can  be 
followed  using  the  increased  fluorescence  of  the  2-aminopurine  probe  which  is  adjacent  to  the  uracil  and  is  very  sensitive  to  base  stacking 
(see  Figure  2)  (5,  6).  Formation  of  P  is  discretely  monitored  by  following  the  decrease  in  tryptophan  fluorescence  of  UDG  that  accompanies 
base  flipping  (5,  6). 


rotation,  a  bent  intermediate  may  form  (lower  right  comer. 
Figure  IB)  before  the  hydrogen  bonds  to  the  target  base  are 
broken  and  the  base  is  expelled  into  the  enzyme  active  site 
(upper  right  comer,  Figure  IB).  A  coupled  reaction  coordi¬ 
nate  for  base  flipping  suggests  that  if  discrete  alterations  in 
the  structural  or  dynamic  features  of  the  DNA  substrate  are 
made,  then  progress  along  one  or  both  coordinates  could  be 
perturbed  in  a  systematic  way.  This  would  allow  a  linear 
free  energy  perturbation  (LFEP)1  analysis  analogous  to  that 
employed  in  simple  chemical  reactions  involving  coupled 
processes  such  as  bond  formation  to  a  nucleophile  and  bond 
breakage  to  a  leaving  group  (13). 

A  LFEP  approach  that  we  are  continuing  to  explore  for 
the  DNA  repair  enzyme  uracil  DNA  glycosylase  (UDG)  is 
to  use  DNA  substrates  in  which  the  number  of  hydrogen 
bonds  between  the  target  uracil  (U)  and  its  opposing  base 
(X)  are  systematically  varied  (Figure  2)  (14).  Our  original 
hypothesis  was  that  removal  of  base  pair  hydrogen  bonds 
would  thermodynamically  facilitate  uracil  flipping  and 
enzyme  binding  by  destabilizing  the  uracil  in  the  DNA  base 
stack  (i.e.,  by  raising  the  free  energy  of  the  lower  left  comer 
of  the  reaction  coordinate  diagram  in  Figure  IB).  A  LFEP 
analysis  has  confirmed  this  initial  expectation,  where  a  strong 
linear  correlation  was  found  between  the  free  energy  of  UDG 
binding  to  a  series  of  DNA  duplexes  with  increasingly 
destabilized  U'X  base  pairs  (m  =  —0.37,  R 2  =  0.86)  (14). 
A  key  question  arising  from  this  thermodynamic  study  is 
the  mechanistic  origin  of  the  dramatically  increased  binding 
affinity  to  destabilized  damaged  sites. 

Here  we  use  a  similar  LFEP  approach  to  correlate  the 
effects  of  stepwise  ablation  of  hydrogen  bonds  in  the  U-X 


1  Abbreviations:  LFEP,  linear  free  energy  perturbation;  2-AP, 
2-aminopurine;  FRET,  fluorescence  resonance  energy  transfer;  M, 
4-methylindole  nucleotide;  D,  diaminopurine  nucleotide;  Y,  pyrene 
nucleotide;  N,  purine  nucleotide;  <I>,  tetrahydrofuran  abasic  nucleotide; 
UF,  2'-/3-fluoro-2'-deoxyuridine  nucleotide;  UDG,  uracil  DNA  glyco¬ 
sylase. 
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Figure  2:  Sequences  of  1 9-mer  DNA  substrates  used  in  the  kinetic 
studies  and  hydrogen-bonding  structures  of  the  UF-X  base  pairs 
based  on  the  literature:  UF-D  (57,  58),  UF-G  (59),  UF-N  (47,  48), 
and  UF,M  (29).  Abbreviations:  M,  4-methylindole  nucleotide;  D, 
diaminopurine  nucleotide;  N,  purine  nucleotide;  UF,  2'-/?-fluoro- 
2'-deoxyuridine  nucleotide. 


base  pair  with  the  rate  of  base  extrusion  (nucleotide  rotation 
reaction  coordinate)  and  the  rate  of  conformational  changes 
in  the  DNA  (bending  reaction  coordinate).  To  monitor  the 
progress  along  these  coordinates,  we  employ  two  indirect 
signal  changes.  The  first  follows  the  increase  in  fluorescence 
of  a  2-aminopurine  (2-AP)  nucleotide  adjacent  to  the  target 
uracil  (5).  This  signal  is  sensitive  to  the  extent  of  stacking 
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of  2-AP  with  the  uracil  (15,  16)  and  thus  reports  on 
conformational  changes  in  the  DNA  as  well  as  progress  along 
the  base  rotation  coordinate  (6,  17—19).  The  second  signal 
involves  a  decrease  in  UDG  tryptophan  fluorescence  that 
accompanies  a  conformational  change  in  UDG  as  it  closes 
around  the  fully  extrahelical  uracil.  This  signal  change  has 
been  shown  by  structural  (5),  kinetic  (5,  6,  17,  20),  and 
mutagenesis  studies  (6,  17)  to  correlate  exclusively  with 
formation  of  the  final  hydrogen-bonding  and  stacking 
interactions  of  the  uracil  base  within  the  enzyme  active  site 
(i.e.,  formation  of  the  final  extrahelical  state  depicted  in  the 
upper  right  comer  of  Figure  IB).  These  studies  suggest  a 
base  flipping  pathway  in  which  DNA  bending  precedes  base 
rotation,  thereby  opening  an  unhindered  passage  by  which 
uracil  may  exit  the  duplex.  Surprisingly,  binding  studies  of 
both  rigid  and  flexible  DNA  duplexes  indicate  that  a  large 
amount  of  enzyme  binding  energy  is  consumed  during  the 
process  of  DNA  bending.  This  finding  appears  to  be  a  general 
feature  of  base  flipping  enzymes  and  is  consistent  with 
previous  suggestions  based  on  structural  and  biochemical 
observations  with  UDG  as  well  as  other  base  flipping 
enzymes  (2,  21—27).  Accordingly,  strong  binding  of  UDG 
to  flexible  target  sites  results  not  only  from  weakened 
hydrogen  bonding  of  the  uracil  but  also  because  flexible  sites 
require  less  enzyme  binding  energy  to  bend.  These  unique 
features  of  target  site  recognition  by  UDG  may  be  shared 
by  other  base  flipping  enzymes  (28). 

EXPERIMENTAL  PROCEDURES 

Materials.  The  2'-deoxynucleoside  phosphoramidites,  CPG 
supports,  and  DNA  synthesis  reagents  were  purchased  from 
Glen  Research  (Sterling,  VA),  except  for  2'-/3-fluoro-2'- 
deoxyuridine  (UF),  4-methylindole  nucleoside  (M),  and 
pyrene  nucleoside  (Y),  which  were  synthesized  as  described 
previously  (5, 29, 30).  The  oligonucleotides  were  synthesized 
using  standard  phosphoramidite  chemistry  on  an  Applied 
Biosystems  392  synthesizer.  The  oligonucleotides  were 
purified  by  anion-exchange  HPLC  (Zorbax  Oligo),  followed 
by  C-18  reversed-phase  HPLC  (Phenomenex  Aqua  column). 
Fractions  containing  pure  oligonucleotide  were  concentrated 
to  dryness  in  vacuo,  redissolved  in  MilliQ  water,  and  stored 
at  —20  °C.  The  purity  of  the  oligonucleotides  was  assessed 
by  matrix-assisted  laser  desorption  mass  spectroscopy  and 
denaturing  polyacrylamide  gel  electrophoresis.  The  concen¬ 
tration  of  each  oligonucleotide  was  determined  using  its 
extinction  coefficient  at  260  nm  (31).  DNA  duplexes  were 
hybridized  in  10  mM  Tris-HCl  (pH  8.0)  and  60  mM  NaCl 
as  described  previously  (5).  The  purification  of  Escherichia 
coli  UDG  has  been  described  previously  (32). 

Kd  Measurements.  The  KD  values  for  binding  of  the 
UF,X  duplexes  to  UDG  were  measured  using  a  kinetic 
competitive  inhibition  fluorescence  assay  under  conditions 
where  the  apparent  K-,  is  equal  to  the  KD  value  (i.e.,  [S]  « 
Km,  where  S  exhibits  rapid  equilibrium  binding).  Reaction 
mixtures  (148.5  pL)  containing  10  mM  Tris-HCl  (pH  8.0), 
60  mM  NaCl,  12.5  pg/mL  BSA,  2  pM  AUPA  substrate,  and 
a  variable  amount  of  the  UF*X  duplex  were  placed  in  a  0.3 
cm  quartz  cuvette,  and  its  fluorescence  emission  at  370  nm 
was  monitored  every  10  s  at  25  °C  on  a  SPEX  FluoroMax-3 
fluorometer  (Acx  =  320  nm)  until  the  signal  stabilized.  UDG 
(1.5  pL)  was  then  added  to  the  reaction  mixture  to  give  a 
final  concentration  of  0.25  nM.  The  progress  of  the  reaction 
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was  then  monitored  as  described  above  for  5  min.  The  Kd 
for  each  duplex  was  determined  by  fitting  to  the  equation 

kJk0  —  1/(1  +  [UF-X]/ATD)  (1) 

where  k\  is  the  inhibited  rate  and  k0  is  the  rate  in  the  absence 
of  competitor  DNA. 

Stopped-Flow  Fluorescence.  The  observed  rate  constants 
for  the  formation  of  UDG-DNA  complexes  were  measured 
on  an  Applied  Photophysics  720  stopped-flow  fluorometer 
(Surrey,  U.K.).  All  measurements  were  performed  under 
pseudo-first-order  conditions  where  the  concentration  of  the 
unlabeled  component  was  at  least  4-fold  greater  the  con¬ 
centration  of  the  labeled  species.  All  measurements  were 
made  using  a  buffer  containing  10  mM  Tris-HCl  (pH  8.0), 
60  mM  NaCl,  and  1  mM  DTT.  For  experiments  where 
changes  in  2-AP  fluorescence  were  observed,  excitation  was 
at  3 15  nm,  and  a  360  nm  long-pass  emission  filter  was  used. 
In  experiments  where  tryptophan  fluorescence  was  measured, 
an  excitation  wavelength  of  290  nm  was  used  with  a  335 
nm  long-pass  emission  filter.  All  kinetic  traces  were  well 
fitted  to  a  first-order  rate  expression  to  obtain  the  observed 
rate  constant  (&0bsd,  eq  2). 

F,  =  AF  exp(  1  -  kobsdt)  +  F0  (2) 

For  the  experiments  in  which  the  rate  of  2-AP  signal 
increase  of  the  DNA  was  monitored,  a  solution  containing  a 
variable  amount  of  UDG  was  rapidly  mixed  with  a  fixed 
concentration  of  2-AP-labeled  DNA.  Final  concentrations  of 
enzyme  and  DNA  after  mixing  were  in  the  range  0.4-6.4 
/rM  and  100  nM,  respectively.  For  the  experiments  in  which 
the  rate  of  tryptophan  signal  decrease  of  the  enzyme  was 
monitored,  a  solution  containing  a  variable  amount  of  DNA 
was  rapidly  mixed  with  a  fixed  concentration  of  the  enzyme. 
In  these  experiments,  the  final  concentrations  of  DNA  and 
enzyme  after  mixing  were  in  the  range  0.4—12.8  «M  and 
1 00  nM,  respectively.  Measurements  with  the  UF,M  duplex 
using  2-AP  fluorescence  could  not  be  accurately  made  above 
1 .6  pM  because  only  20%  of  the  observable  signal  remained 
after  the  ~1  ms  dead  time  of  the  instrument  had  elapsed. 
As  previously  observed,  plots  of  kobsd  against  [UDG]  or 
[UF-X]  were  hyperbolic,  indicating  a  multistep  binding 
mechanism.  The  kinetic  parameters  were  extracted  by  fitting 
the  data  to  a  minimal  two-step  binding  model  (eqs  3  and  4), 

E  +  S  =f=^  ES  sf=S:  E*F  (3) 

k~  |  k~2 

= ys]  +  *„fr  (4) 

obsd  K'[  S]  +  1  W 

where  K'  —  k\/(k-\  +  kmax)  is  the  apparent  affinity  constant, 
kon  —  K'(ki  +  k~i)  is  the  apparent  second-order  association 
rate  constant,  and  k0^=  k-\k-2/(k-\  +  k2  +  k-i)  is  the  overall 
dissociation  rate  constant.  The  maximum  rate  constant  for 
the  unimolecular  rearrangements  detected  by  the  2-AP  and 
tryptophan  fluorescence  signals  is  measured  by  the  asymptote 
(or  £maxT,p)  =  k2  +  k-2  (5).  Although  base  flipping 
has  been  previously  shown  to  involve  two  internal  steps 
rather  than  the  one  shown  in  eq  3  (6, 1 7, 20),  the  simplified 
analytical  expression  of  eq  4  is  very  useful  for  comparing 
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the  kinetic  behavior  of  a  series  of  substrates  or  mutated 
enzymes  (6). 

Very  slow  dissociation  rate  constants  determined  from  the 
concentration  dependence  of  ft0bsd  (eq  3)  have  a  high  degree 
of  uncertainty  because  fr0fr  is  derived  from  the  v-intercept  of 
a  plot  of  kobsi  against  [S].  Under  such  conditions,  a  trapping 
experiment  is  required  to  obtain  a  more  precise  value  for 
k0ff.  Trapping  experiments  for  the  substrates  containing  D, 
A,  G,  and  N  in  the  base  pair  with  uracil  were  performed  by 
rapidly  mixing  a  solution  consisting  of  UDG  (200—500  nM) 
and  2-AP-labeled  DNA  (200  nM)  with  a  large  molar  excess 
of  nonfluorescent  trapping  duplex  (20  uM  AUFA-11)  (6) 
using  a  stopped-flow  fluorometer.  The  time-dependent 
decrease  in  2-AP  fluorescence  as  the  bound  DNA  irreversibly 
dissociated  was  fit  to  a  single-exponential  decay  to  obtain 
k0f{  (eq  5).  Because  of  the  very  slow  dissociation  of  UF-M, 

F,  —  AF  exp  (-kof(t)  +  F0  (5) 

its  ftoff  was  measured  by  manually  mixing  a  solution 
consisting  of  UDG  (100  nM)  and  PUF,M  (100  nM)  with  1 
yM  high-affinity  trap,  AUFA/TMT-15  (14).  This  reaction 
was  performed  using  a  0.3  cm  path  length  quartz  cuvette 
(150  «L)  and  a  Spex  FluoroMax-3  fluorometer.  The  time- 
dependent  decrease  in  2-AP  emission  at  370  nm  was 
followed  with  excitation  at  3 1 5  nm. 

RESULTS 

Binding  of  UDG  to  Destabilized  Damaged  Sites.  We 
synthesized  a  series  of  19-mer  DNA  duplexes  containing  a 
series  of  UF-X  base  pairs  using  a  sequence  based  on  a  duplex 
previously  used  in  rapid  kinetic  studies  of  base  flipping  by 
UDG  (Figure  2)  (5).  The  2'  fluorinated  deoxyuridine 
substrate  analogue  (UF)  is  an  extremely  slow  substrate  for 
UDG  (/ 1/2  ~  1  day),  allowing  measurements  of  DNA  binding 
and  base  flipping  without  the  complication  of  glycosidic  bond 
cleavage  (5).  The  affinity  of  UDG  for  these  duplexes  was 
measured  using  a  competitive  inhibition  kinetic  assay  (14). 
As  previously  observed  for  a  series  of  15-mer  duplexes  with 
a  different  sequence  (14),  the  Kq  values  for  the  19-mers 
decreased  incrementally  as  the  UF*X  base  pair  weakened 
(Figure  3A  and  Table  1).  The  UF-M  duplex,  which  contains 
no  hydrogen  bonds,  binds  a  striking  ~43400-fold  (6.3 
kcal/mol)  more  strongly  than  the  duplex  containing  the 
UF*D  base  pair  with  three  hydrogen  bonds.  Of  note,  the 
incremental  decrease  in  the  free  energy  of  binding  was  only 
about  1  kcal/mol  when  the  first  and  second  hydrogen  bonds 
were  removed  from  the  UF>D  base  pair  (Figure  3B),  which 
is  similar  to  estimates  of  the  free  energy  contribution  of 
individual  hydrogen  bonds  to  the  stability  of  duplex  nucleic 
acids  (53).  In  contrast,  removal  of  the  last  hydrogen  bond  to 
form  the  UF,M  pair  resulted  in  a  much  larger  4  kcal/mol 
decrease  in  binding  free  energy,  indicating  the  presence  of 
additional  energetic  contributions.  These  additional  contribu¬ 
tions  may  involve  a  loss  of  both  hydrogen  bond  and  base 
stacking  interactions,  leading  to  a  large  increase  in  confor¬ 
mational  entropy  of  the  UF-M  base  pair  (14).  This  result 
suggests  that,  with  respect  to  the  free  energy  of  binding, 
increased  flexibility  of  the  site  may  play  a  more  important 
role  than  the  enthalpic  benefit  of  removing  single  hydrogen 
bonds  (see  below). 


H-bonds 


Figure  3 

Figure  3:  Binding  affinities  of  UDG  for  UF-X  DNA  duplexes  and 
incremental  change  in  binding  free  energy  when  each  hydrogen 
bond  in  removed.  (A)  Trend  in  dissociation  binding  constants  as 
base  pair  hydrogen  bonds  are  removed  (see  also  Table  1).  The 
number  of  hydrogen  bonds  in  the  UF-X  base  pair  is  shown  to  the 
right  of  the  bars.  (B)  Incremental  change  in  binding  free  energy  as 
each  hydrogen  bond  is  removed  from  the  UF,X  base  pair  (AAGbind 
=  — ftrin  Kd"/Ko"~\  where  n  is  the  number  of  hydrogen  bonds  in 
the  base  pair). 

Binding  Kinetics.  To  dissect  the  origins  of  the  dramatically 
enhanced  binding  affinity  of  the  duplexes  with  destabilized 
UF,X  base  pairs,  the  kinetics  of  association  and  dissociation 
were  measured.  Association  rates  were  measured  using 
stopped-flow  fluorescence  measurements  by  monitoring 
either  the  increase  in  2-AP  fluorescence  that  accompanies 
uracil  unstacking  and  DNA  bending  (Figure  4A)  or  the 
decrease  in  tryptophan  fluorescence  that  marks  the  attainment 
of  the  final  extrahelical  state  (Figure  4B)  (5,  6,  17).  As 
observed  in  previous  studies  using  duplexes  with  UF,A  and 
UF-G  base  pairs  (5),  plots  of  k0bsi  versus  concentration  were 
hyperbolic  when  either  the  2-AP  or  tryptophan  fluorescence 
signals  were  monitored  (panels  A  and  B  of  Figure  5, 
respectively).  This  kinetic  behavior  indicates  a  change  in  rate- 
limiting  step  from  bimolecular  encounter  at  low  concentra¬ 
tions  of  the  varied  reactant  to  a  unimolecular  conformational 
change  of  the  DNA  and  enzyme  at  high  reactant  concentra¬ 
tions.  All  hyperbolic  plots  were  fit  to  the  two-step  binding 
model  (eq  4),  and  the  kinetic  constants  obtained  from  this 
analysis  are  reported  in  Table  1. 

There  are  revealing  aspects  of  the  kinetics  for  the  2-AP 
and  tryptophan  fluorescence  changes.  First,  for  the  substrate 
with  the  most  stable  base  pair  (UF,D),  the  association  kinetics 
measured  using  the  2-AP  signal  (Zr0n2AP)  is  3-fold  faster  than 
when  the  tryptophan  signal  is  followed  (konTrp)  (Table  1). 


Extrahelical  Base  Recognition 
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Table  1 :  Binding  Affinities  and  Kinetic  Parameters  for  UDG  Association  and  Dissociation  with  UF-X  Duplexes 

b  2AP  b  Trp 

Aon  Aon 


duplex 

XD(nM) 

K'2M 

(uM-1) 

K'  Trp(«M-1) 

-1. 

s'1) 

(«M 

-1. 

s1) 

k„nis  ’) 

XDcalc(nM)‘' 

UF-D 

434  ±  100“ 

0.29 

± 

0.05 

0.12 

± 

0.02 

206 

± 

22 

72 

± 

7 

44  ±  1.3 

214 

± 

24 

UF-A 

80  ±  20"'’ 

0.08 

± 

0.01 

0.13 

± 

0.01 

236 

± 

10 

172 

± 

8 

26  ±  0.3 

110 

± 

4 

UF-G 

50  ±  20"'’ 

0.10 

± 

0.01 

0.25 

± 

0.08 

292 

± 

13 

331 

± 

56 

6.9  ±  0.04 

24 

± 

1 

UF-N 

8.2  ±  0.3" 

0.20 

± 

0.04 

0.14 

± 

0.01 

282 

± 

28 

315 

± 

5 

6.0  ±  0.05 

21 

± 

2" 

UF-M 

0.01  ±0.002" 

0.32 

± 

0.15 

0.64 

± 

0.08 

292 

± 

40 

666 

± 

45 

0.018  ±0.001 

0 

062  ±  0.009 

°  Ko  was  determined  by  directly  monitoring  binding  of  DNA  to  UDG  using  2-AP  fluorescence  (see  Experimental  Procedures). b  Values  previously 
reported  (5). "  K0  was  determined  using  competitive  inhibition  assay  with  AUF-X  substrate  (see  Experimental  Procedures).  The  affinity  of  AUF-A 
is  2-fold  greater  than  that  of  PUF-A  (17).  J  Ka  calculated  from  Ar0fAn2AP. 


A 


Figure  4:  Stopped-flow  fluorescence  kinetic  measurements  of 
UDG  association  with  UF,X  duplex  DNA.  (A)  Approach  to 
equilibrium  association  rate  of  PUF-N  (0. 1  /rM)  and  UDG  (0.4  ,uM) 
monitored  by  the  increase  in  2-AP  fluorescence.  (B)  Approach  to 
equilibrium  association  rate  of  UDG  (0.1  /rM)  and  AUF,N  (0.4  /<M) 
monitored  by  the  decrease  in  tryptophan  fluorescence  of  UDG. 

However,  as  hydrogen  bonds  are  removed  from  the  base  pair 
kon2AP  and  k0 nTrp  become  indistinguishable.  This  result  arises 
because  k0 n2AP  is  nearly  invariant  across  the  whole  series  of 
duplexes,  while  k0 nTrp  increases  by  about  9-fold  (Figure  5C).2 
The  observation  that  klmTrp  is  slower  than  kon2M  when  the 
base  pair  is  strong  requires  that  the  conformational  change 
in  the  enzyme  lags  behind  the  structural  perturbation  in  the 
DNA  that  results  in  the  2-AP  fluorescence  increase.  The 
strong  dependence  of  konj!v  on  base  pair  hydrogen  bonding 
indicates  that  these  hydrogen  bonds  are  broken  very  late 


2  The  asymptotic  k0 bsd  values  for  the  2-AP  signal  (£ma*2AP  =  fo2AP  + 
k- 22AP;  see  eq  3)  increase  by  about  2-fold  between  UF-D  and  the  UF,A 
or  UF-G  substrates  but  then  decrease  from  these  peak  values  by  about 
2-fold  for  the  UF-N  and  UF-M  substrates.  This  behavior  may  be 
qualitatively  understood  using  the  simplified  two-step  kinetic  analysis 
we  have  employed.  That  is,  £22AP  and  £-22AP  are  complex  rate  constants 
that  reflect  the  interconversion  of  Ij,  E,  and  P  on  the  enzyme:  £,m„2AP 
first  increases  as  the  conversion  of  I2  — *  P  becomes  more  rapid  due  to 
hydrogen  bond  ablation  and  then  decreases  because  the  reverse  rate 
(i_22AP)  becomes  negligible  for  the  UF-N  and  UF-M  substrates  due  to 
stabilization  of  I2  and  P  (i.e.,  kmm2A?  ~  &22AP;  see  Discussion). 


♦  PUf/D 
a  puf/A 
■  PUf/G 
▼  PUf/N 

•  PUf/M 


♦  auf/d 

A  AUF/A 
■  AUF/G 
▼  AUF/N 

•  AUF/M 


2-aminopurine  tryptophan 

fluorescence  fluorescence 

Figure  5:  Concentration  dependence  of  the  apparent  association 
rate  constants  (Arobs)  for  the  five  different  UF,X  base  pairs  and 
relative  kon  values  as  compared  to  the  UF-D  substrate.  (A) 
Dependence  of  kobs  on  the  concentration  of  UDG  followed  by 
changes  in  2-AP  fluorescence.  Data  were  fit  to  the  two-step  kinetic 
model  (eqs  3  and  4).  (B)  Dependence  of  /robs  on  the  concentration 
of  DNA  followed  by  changes  in  Trp  fluorescence.  Data  were  fit  to 
the  two-step  kinetic  model  (eqs  3  and  4).  (C)  Effect  of  changing 
UF-X  base  pair  strength  on  k0T],  expressed  as  the  ratio  of  Aron(UF-X) 
to  £on(UF,D).  Ratios  determined  from  the  kinetics  of  the  2-AP  and 
tryptophan  fluorescence  changes  are  displayed  separately. 


along  the  reaction  coordinate,  after  the  step  that  is  reported 
on  by  the  2-AP  signal  and  during  the  final  conformational 
step  when  the  enzyme  closes  around  the  uracil  base.  This 
result  implies  that  the  6-fold  2-AP  signal  increase  reflects 
the  formation  of  an  intermediate  (I2,  Figure  1)  in  which  the 
uracil  stacking  interactions  with  the  adjacent  2-AP  base  are 
perturbed,  but  the  uracil  retains  hydrogen  bonding  with  its 
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Figure  6:  Determination  of  k0rr  using  a  trapping  experiment.  A 
complex  of  PUF-N  (0.20  fiM)  and  UDG  (0.25  /(M)  was  rapidly 
mixed  with  a  large  excess  of  nonfluorescent  AUFA-I1  single- 
stranded  DNA  (20  /rM)  (6).  The  k0r;  is  determined  from  fitting  the 
time-dependent  decrease  in  the  2-AP  fluorescence  of  PUF*N  as  it 
is  irreversibly  released  from  UDG  (eq  5). 

partner  base  (6,  17).  Since  the  preceding  intermediate  (Ii, 
Figure  1)  is  not  detectably  bent  (7),  and  the  subsequent 
extrahelical  state  is  bent  by  about  40°  (8),  then  the  structural 
implication  of  the  increase  in  2-AP  fluorescence  is  that  the 
DNA  is  bent  in  the  transition  state  leading  to  /2 ,3  although 
the  magnitude  of  the  bending  at  this  point  along  the  reaction 
coordinate  is  not  known.  It  is  important  to  point  out  that  the 
9-fold  increase  in  kon  pales  in  comparison  with  the  ap¬ 
proximately  43000-fold  increase  in  UDG  binding  affinity 
when  these  three  hydrogen  bonds  are  ablated.  Thus,  the  vast 
majority  of  the  binding  affinity  increase  is  not  attributable 
to  kon. 

The  above  results  strongly  suggest  that  the  large  thermo¬ 
dynamic  benefit  of  removing  hydrogen  bonds  in  the  UF-X 
base  pair  arises  from  profound  differences  in  the  dissociation 
rates  of  these  duplexes.  To  establish  this  point,  the  dissocia¬ 
tion  rates  (/:0nj  were  measured  using  an  irreversible  trapping 
experiment  (5)  in  which  a  solution  of  UDG-(PUF,X)  was 
rapidly  mixed  with  a  large  excess  of  unlabeled  trapping 
DNA,  and  the  decrease  in  the  2-AP  signal  was  followed 
(Figure  6).  The  dissociation  rate  constants  are  reported  in 
Table  1  and  are  consistent  with  the  y-intercepts  (Aofy)  esti¬ 
mated  from  the  concentration  dependence  of  the  association 
rates  in  Figure  5A,B.  The  dissociation  rates  were  strongly 
dependent  on  the  nature  of  the  UF-X  base  pair,  decreasing 
by  over  4  orders  of  magnitude  over  the  entire  series,  with 
A:off(UF-D)  =  44  s">  and  MUF-M)  =  0.018  s“l.  These  data 
provide  strong  evidence  that  removal  of  target  base  pair 
hydrogen  bonds  leads  to  profound  stabilization  of  one  or 
more  enzyme-bound  species  on  the  uracil  flipping  pathway. 
Since  the  kinetic  effect  of  removing  hydrogen  bonds  appears 
very  late  during  the  uracil  flipping  process  (i.e.,  during  the 


3  We  were  unable  to  detect  DNA  bending  by  UDG  employing 
standard  fluorescence  resonance  energy  transfer  (FRET)  methods  and 
DNA  substrates  with  fluorescence  donor  and  acceptor  groups  on  each 
end  (51,  52).  We  attribute  this  result  to  interaction  of  the  rhodamine 
fluorophore  with  binding  sites  in  the  free  DNA  (53,  54),  which  may 
be  disrupted  when  the  enzyme  binds,  resulting  in  no  FRET  change. 
Based  on  crystallographic  models  of  bent  DNA  bound  to  UDG  (not 
shown),  we  would  have  expected  a  substantial  1.4-fold  FRET  increase. 
The  same  negative  result  was  obtained  when  the  donor  or  acceptor 
was  placed  on  the  3'  or  5'  end  of  the  uracil-containing  strand  or  when 
a  15-mer  or  19-mer  duplex  was  used.  We  note  that  an  unexpected  and 
small  15%  FRET  decrease  was  recently  reported  for  £coRI  methyl- 
transferase  which  also  flips  a  DNA  base  by  a  putative  bending 
mechanism,  although  an  increase  was  expected  (35).  We  conclude  that 
FRET  results  should  be  interpreted  cautiously  in  base  flipping  systems. 


P  Kd 


Figure  7:  Logarithmic  dependences  of  km  and  Xr0fr  on  the  binding 
affinity  of  UDG  for  UF-X  DNA  (pATo).  (A)  Logarithmic  dependence 
of  &on  on  pXD.  The  data  were  fit  to  the  empirical  equation  log  k0 „ 
=  log  Liifr  —  Iog(  1  +  1  Op^-'P^d),  where  is  the  rate  constant  for 
diffusional  encounter  and  pX0  is  an  empirical  constant.  (B)  Linear 
dependence  of  log  k0n  on  binding  affinity  (slope  =  0.74  ±  0.07, 
R 2  =  0.97).  For  reference,  the  number  of  hydrogen  bonds  in  each 
destabilized  UF,X  base  pair  is  indicated  above  each  data  point. 

last  conformational  step  monitored  by  tryptophan  fluores¬ 
cence;  see  above),  then  the  stabilized  species  must  be  one 
or  both  of  the  complexes  shown  in  the  upper  and  lower  right- 
hand  comers  of  Figure  IB  (I2,  P). 

Linear  Free  Energy  Correlations.  The  quantitative  de¬ 
pendence  of  kon  and  kofr  on  UDG  binding  affinity  provides 
further  insights  into  the  mechanistic  basis  for  the  enhanced 
binding  of  duplexes  with  destabilized  base  pairs.  A  plot  of 
log  k0„2AP  versus  pXD  is  nearly  flat  across  the  whole  series 
(Figure  7A,  triangles)  with  an  average  value  kon  —  2.8  x 
108  M_1  s~'  indicative  of  near  diffusion-controlled  binding 
except  for  the  substrate  with  the  most  stable  base  pair 
(UF-D)  which  falls  slightly  below  this  average  value.  Thus, 
all  internal  kinetic  steps  preceding  and  including  the  2-AP 
fluorescence  change  are  rapid  as  compared  to  the  encounter 
rate.  In  contrast,  k0 nTrp  shows  a  linear  increase  as  a  function 
of  p Kd  for  the  more  stable  UF-D  and  UF-A  duplexes, 
followed  by  downward  curvature  to  a  plateau  value  of  ft0nTip 
=  4  x  108  M-1  s“‘  for  the  most  destabilized  base  pairs 
(Figure  7A,  squares).  This  behavior  indicates  a  change  in 
rate-limiting  step  as  a  function  of  base  pair  stability  from 
the  step  giving  rise  to  the  tryptophan  fluorescence  change 
to  that  involving  diffusion-controlled  encounter.  We  conclude 
that  the  thermodynamic  stability  of  the  UF-X  base  selectively 
alters  the  activation  barrier  for  the  step  reported  by  k0 nTrp. 

In  contrast  with  log  kon,  there  is  a  linear  dependence  of 
log  k0ff  on  pAT>,  with  the  dissociation  rate  decreasing  steeply 
as  UDG  affinity  increases  (Figure  7B).  From  the  slope  of 
this  correlation  (—0.74  ±  0.07),  it  can  be  inferred  that  about 
74%  of  the  difference  in  binding  energies  between  any  two 
UF*X  duplexes  arises  from  k0ff.  The  simplest  mechanism  that 
accounts  for  the  small  effect  of  base  pair  strength  on  kon, 
and  the  large  effect  on  k0 fr,  is  the  progressive  stabilization 
of  one  or  more  enzyme-bound  species  as  hydrogen  bonds 
are  removed  (see  Discussion). 
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Table  2:  Effect  of  Pyrene  Nucleoside  on  Binding  Affinity  of  UDG 
to  Uracil  and  the  Abasic  Site  Containing  DNA" 

Duplex  DNA  KD  (nM) 


u  F^f° 

N 

DNA'  J 

-H-N  u 

300  ±  50 

UF 

v, 

120  ±30 

uW° 

L  M  H 

°NA'fT 

"Ip: 

0.24  ±0.03 

4> 

DNA' 

H 

H-N  N 

ffi*.  A 

\N  DNA 

100  ±10 

a> 

DNA 

eft? 

8700  ±  675 

DNA 

DNA 

H3C  _ 

yp,  m 

W/  DNA 

18  ±3.0 

"The  sequence  of  DNA  corresponds  to  15-mer  duplexes  in  ref  14. 
<t>  =  tetrahydrofuran  abasic  nucleotide  analogue. 


Base  Pair  Hydrogen  Bonds,  DNA  Flexibility,  and  UDG 
Binding.  To  test  the  role  of  DNA  flexibility  on  binding 
affinity,  we  employed  a  novel  structure-reactivity  approach. 
This  strategy  takes  advantage  of  the  decrease  in  flexibility 
of  the  target  site  when  it  is  modified  with  a  rigid  pyrene 
nucleotide  on  one  strand  (Table  2)  (17,  30,  34-36).  We 
originally  used  a  uracil-pyrene  base  pair  (UF*Y)  to  preor¬ 
ganize  uracil  in  an  extrahelical  conformation  (34),  with  the 
rationale  that  the  planar  aromatic  ring  structure  of  pyrene 
would  fully  encompass  the  entire  volume  normally  occupied 
by  the  standard  UF-A  base  pair,  thereby  presenting  the  uracil 
to  the  enzyme  in  an  extrahelical  conformation  that  favors 
binding  (i.e.,  the  energetic  penalty  for  uracil  flipping  is 
prepaid,  allowing  for  tighter  binding).  In  accord  with  this 
reasoning,  DNA  that  contained  a  UF,Y  pair  fully  rescued 
the  catalytic  and  base  flipping  activities  of  UDG  mutants 
that  were  incapable  of  pushing  uracil  from  the  DNA  base 
stack  in  the  context  of  a  UF-A  base  pair  (17,  34). 

Although  base  flipping  mutants  showed  dramatic  increases 
in  their  site-specific  binding  affinity  and  catalytic  activity 
when  a  UF,Y  pair  was  used,  wild-type  UDG  bound  with 
nearly  equal  affinity  to  DNA  containing  either  a  UF*Y  or 
UF,A  base  pair  (34).  This  finding  is  reaffirmed  in  Table  2 
where  the  binding  affinity  of  UDG  to  a  15-mer  duplex 
containing  a  UF*A  and  UF-Y  pair  is  shown  to  differ  by  only 
2.5-fold.  This  result  is  in  dramatic  contrast  to  the  1200-fold 
tighter  binding  of  the  UF-M  duplex  as  compared  to  UF,A 
using  the  same  15-mer  context  (Table  2).  These  observations 
indicate  that  there  is  a  cryptic  energetic  penalty  for  binding 
of  the  UF-Y  duplex  that  is  not  present  for  binding  of  the 
UF*M  duplex  even  though  UF,M  and  UF,Y  both  lack 
hydrogen  bonds  to  the  uracil  base.  Thus  a  key  question  is 
the  basis  for  the  dramatically  tighter  binding  of  the  UF-M 
duplex  as  compared  to  UF,Y  (A AG  —  3.1  kcal/mol). 

A  likely  explanation  is  that  the  pyrene  nucleoside  serves 
to  stiffen  the  DNA  at  the  apex  of  the  bend  induced  during 
base  flipping  and  thus  requires  a  greater  expenditure  of 


binding  energy  by  UDG  in  order  to  distort  the  DNA  as 
compared  to  UF,M.  This  idea  is  consistent  with  previous 
computational  findings  suggesting  that  decreased  DNA 
flexibility  can  lead  to  reduced  UDG  binding  affinity  by 
increasing  the  resistance  of  the  DNA  to  adopt  the  bend 
enforced  by  protein  binding  (10).  Increased  structural  rigidity 
introduced  by  pyrene  is  chemically  reasonable  because  the 
extended  Jt  aromatic  system  of  the  pyrene  base  pair  mimic 
is  much  more  rigid  than  the  native  UF-A  base  pair  which  is 
only  restrained  by  noncovalent  hydrogen  bonding.  Assuming 
this  interpretation  for  binding  of  the  UF*Y  base  pair,  its 
unfavorable  rigidity  is  offset  by  the  energetic  benefit  of 
extrahelical  preorganization  of  the  uracil.  Thus,  the  binding 
affinity  of  UF*Y  to  UDG  is  similar  to  a  UF-A  base  pair 
because  of  compensatory  energetic  effects,  but  UF,M  binds 
more  tightly  because  of  the  increased  flexibility  of  the  site 
(Table  2). 

To  further  test  the  proposal  that  rigidity  of  the  DNA  has 
an  adverse  effect  on  UDG  binding  affinity,  we  tested  the 
effect  of  the  opposing  base  (A,  M,  or  Y)  on  binding  of  DNA 
that  contained  an  abasic  site  product  analogue  (4>)  (Table 
2).  Structural  studies  have  shown  that  UDG  can  bind 
specifically  to  and  flip  out  abasic  sites  to  form  a  structurally 
indistinguishable  complex  as  compared  to  substrate  analogue 
DNA  (2,  8,  37,  38).  However,  because  there  is  no  require¬ 
ment  for  rupture  of  a  base  pair  upon  flipping  of  an  abasic 
nucleotide,  and  there  is  no  uracil  attached  to  form  hydrogen 
bonds  and  stacking  interactions  in  the  UDG  active  site,  the 
changes  in  binding  affinity  between  the  4>*Y  and  5>‘M 
duplexes  should  largely  reflect  differences  in  the  intrinsic 
flexibility  of  the  target  site.  As  shown  in  Table  2,  4>-M  binds 
~480-fold  more  tightly  than  «J>-Y  and  5-fold  more  tightly 
than  <I>-A.  The  observed  greater  affinity  of  UDG  for  <1>-M 
over  <f>-Y  is  identical  to  that  of  UF-M  and  UF,Y  (Table  2). 
We  conclude  from  these  results  that  pyrene  introduces  an 
unfavorable  free  energy  contribution  to  binding  of  as  much 
as  3.7  kcal/mol  arising  from  its  increased  rigidity  and 
stacking  as  compared  to  sites  that  contain  A  or  M  as  the 
opposing  base.  In  conclusion,  the  observation  that  flexible 
abasic  sites  bind  with  much  higher  affinity  than  rigid  abasic 
sites  provides  strong  experimental  evidence  that  specific 
recognition  involves  a  flexibility  component  regardless  of 
the  presence  of  a  uracil  base. 

DISCUSSION 

The  mechanism  by  which  enzymes  obtain  extraordinary 
specificity  for  their  respective  target  sites  in  a  large  back¬ 
ground  of  random  DNA  sequences  has  long  been  an  active 
area  of  research.  In  this  general  area,  DNA  glycosylases  have 
evolved  a  unique  solution  to  this  problem  that  is  divergent 
from  the  sequence-dependent  mechanisms  of  restriction 
enzymes,  transcription  factors,  and  repressor  proteins.  Ob¬ 
taining  a  fundamental  understanding  of  the  kinetic  and 
thermodynamic  origins  of  DNA  repair  glycosylase  specificity 
has  the  potential  to  allow  rational  targeting  of  these  enzymes 
to  engineered  DNA  sites  that  display  features  that  promote 
enhanced  binding  and  reactivity.  Indeed,  we  have  previously 
shown  that  simple  rules  derived  from  mechanistic  studies 
of  base  flipping  and  catalysis  can  be  used  to  alter  UDG’s 
specificity  to  recognize  a  cytosine  base  opposite  to  pyrene 
(i.e.,  a  OY  base  pair)  rather  than  uracil  (39).  Such  findings 
offer  the  promise  of  altering  the  coding  sequences  of  genes 
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in  vivo  using  engineered  glycosylases  in  combination  with 
assisting  targeting  molecules  such  as  pyrene-containing 
oligonucleotides. 

The  Reaction  Coordinate  for  Base  Flipping.  It  is  difficult 
to  imagine  the  pathway  for  flipping  a  base  180°  from  the 
DNA  base  stack  ultimately  leading  to  its  precise  docking 
into  an  enzyme  active  site.  It  is  equally  difficult  to  view  how 
an  enzyme  could  form  specific  interactions  to  guide  a  target 
base  through  a  trajectory  that  arcs  widely  out  of  the  DNA 
helix  (Figure  1  A)  (7,  55).  The  present  data  suggest  plausible 
mechanistic  answers  to  these  questions  that  are  consistent 
with  previous  NMR  dynamic  and  structural  studies  of  the 
earliest  steps  in  the  base  flipping  process  (7,  55)  and  also 
crystallographic  studies  of  the  final  extrahelical  complex  (8). 

To  understand  how  an  enzyme  might  facilitate  the  overall 
base  flipping  process,  it  is  useful  to  consider  the  intrinsic 
kinetic  and  thermodynamic  problems  that  must  be  overcome 
in  extruding  a  base  from  the  B  DNA  base  stack  (Figure  8). 
First,  NMR  imino  proton  exchange  experiments  have  estab¬ 
lished  that  T*A,  T*G,  and  U-A  base  pairs  rapidly  open  at 
room  temperature  (7,  40,  41).  These  opening  rates  are  greater 
than  or  equal  to  the  rate  constants  for  kinetic  steps  on  the 
base  flipping  of  UDG  (Table  1),  indicating  that  spontaneous 
base  pair  opening  provides  kinetically  competent  motions 
for  seeding  the  enzymatic  base  flipping  process.  Indeed, 
recent  studies  of  imino  proton  exchange  of  UDG-bound  DNA 
have  established  that  UDG  does  not  alter  the  opening  rate 
but  instead  substantially  slows  the  closing  rate  constant  of 
open  T*A  base  pairs  (7).  This  result  is  intuitively  satisfying 
because  the  base  pair  closing  rates  in  B  DNA  are  exceedingly 
fast  (~107  s"1)  (42, 43),  and  one  major  problem  in  facilitating 
the  migration  of  a  base  forward  along  the  flipping  trajectory 
is  to  prevent  its  retrograde  motion  back  into  the  DNA  base 
stack  (Figure  8A).  On  the  basis  that  adjacent  T-A  and  G'C 
base  pairs  showed  identical  closing  rates  in  the  UDG 
complex  (MO5  s-1),  but  these  rates  differed  by  25-fold  in 
the  free  DNA,  we  proposed  that  UDG  uses  nonspecific  DNA 
backbone  interactions  to  increase  the  lifetime  of  a  high- 
energy  open  state  of  the  bound  DNA  that  occurs  very  early 
on  the  flipping  trajectory  (Ii,  Figure  8A)  (7).  This  open  state 
of  UDG-bound  DNA  is  still  extremely  unstable  and  exists 
in  dynamic  equilibrium  with  closed  B  DNA.  Nevertheless, 
the  equilibrium  is  now  shifted  more  toward  the  open  state 
as  compared  to  the  free  DNA  (Figure  8A). 

The  next  intermediate  (I2)  on  the  flipping  pathway  has  an 
altered  B  DNA  structure  such  that  the  uracil  stacking 
interactions  with  the  adjacent  2-AP  are  perturbed,  but  on 
average,  the  base  pair  hydrogen  bonds  are  still  intact  (Figure 
8A).  The  experimental  observations  that  support  these 
features  of  I2  are  (i)  the  rapid  2-AP  fluorescence  change 
indicating  a  perturbation  in  the  stacking  of  uracil  with  2-AP 
during  the  formation  of  I2  (Figure  4A)  and  (ii)  that  the  kinetic 
step  between  I2  and  the  final  extrahelical  state  (P)  is 
kinetically  enhanced  as  hydrogen  bonds  are  removed  from 
the  U-X  base  pair  (Figure  5C),  indicating  that  these  bonds 
persist  after  formation  of  I2.  Once  the  final  state  is  attained 
and  the  uracil  is  docked  and  held  securely  in  the  active  site 
by  stacking  and  hydrogen-bonding  interactions  (8,  38,  44— 
46),  the  extrahelical  base  cannot  reenter  the  DNA  stack 
without  climbing  the  reverse  activation  barrier  back  to  I2. 
Thus,  the  base  flipping  process  may  be  viewed  as  the 
progressive  formation  of  enzyme-bound  states  in  which  the 
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uracil  base  spends  an  ever  increasing  time  in  an  extrahelical 
conformation  as  depicted  in  Figure  8A. 

The  Importance  of  Being  Flexible.  The  kinetic  and 
thermodynamic  data  indicate  that  as  the  base  pair  hydrogen 
bonds  are  removed,  the  second  intermediate  (I2),  or  the  final 
extrahelical  state  (P),  is  increasingly  stabilized  (Figure  8B). 
The  results  also  reveal  that  removal  of  the  first  two  hydrogen 
bonds  alters  binding  affinity  by  M  kcal/mol  per  hydrogen 
bond  but  that  removal  of  the  final  hydrogen  bond  results  in 
a  significantly  larger  4  kcal/mol  increase  in  the  binding 
affinity  (Figure  3B).  Assuming  that  the  base  pair  remains 
largely  intact  when  the  first  two  hydrogen  bonds  are  broken, 
which  is  supported  by  both  NMR  and  potassium  perman¬ 
ganate  oxidation  sensitivity  measurements  (14,  47,  48),  then 
the  1  kcal/mol  free  energy  increment  may  provide  an  estimate 
of  the  individual  hydrogen  bond  enthalpy  (55).  The  free 
energy  change  when  the  last  bond  is  removed  is  substantially 
greater  than  estimates  of  DNA  base  hydrogen  bond  energy 
(55,  49),  indicating  the  loss  of  multiple  energetic  interactions, 
perhaps  including  stacking  of  the  uracil  with  adjacent  bases. 
We  suggest  that  one  ramification  of  total  ablation  of 
hydrogen  bonds  in  the  UF,X  pair  is  increased  conformational 
flexibility  of  the  target  site.  Such  conformational  flexibility 
would  be  expected  to  facilitate  DNA  bending  and,  therefore, 
binding  affinity.  The  conclusion  that  hydrogen  bond  ablation 
increases  target  site  flexibility  is  further  supported  by  the 
3.7  kcal/mol  unfavorable  effect  of  the  rigid  pyrene  (Y) 
nucleotide  on  binding  as  compared  to  4-methylindole  (M) 
when  these  nucleotides  are  placed  opposite  to  uracil  or  an 
abasic  site  in  duplex  DNA  (Table  2).  We  infer  that  pyrene, 
by  forming  strong  stacking  interactions  and  occupying  a 
space  that  spans  the  width  of  the  DNA  duplex,  significantly 
increases  the  resistance  of  the  DNA  toward  bending,  which 
is  in  turn  reflected  in  the  binding  affinity  measurements.  Of 
course,  the  other  destabilized  base  pairs  may  also  exhibit 
changes  in  stacking  energetics,  but  these  differences  are 
expected  to  be  small  given  the  very  conservative  changes  in 
the  nature  of  the  opposing  base  (50).  These  differences  in 
stacking  energies,  if  present,  would  also  be  manifested  in 
the  flexibility  of  the  DNA  and  would  therefore  contribute 
to  the  energetics  of  enzyme-induced  DNA  bending. 

A  free  energy  diagram  depicting  the  significant  effect  of 
target  site  flexibility  on  binding  is  shown  in  Figure  8B.  The 
diagram  highlights  how  weaker  UF-X  base  pairs  are  easier 
to  bend  due  to  their  increased  conformational  flexibility, 
leading  to  a  decrease  in  the  free  energy  of  the  bent  DNA, 
reflected  in  the  stability  of  I2  and  P  in  Figure  8B.4  In  addition 
to  the  energetic  benefit  of  increased  flexibility  of  the  UF-M 
substrate  (~3  kcal/mol),  the  P  state  realizes  the  additional 
energetic  benefit  of  forming  hydrogen  bonds  with  the  uracil 
base  without  the  penalty  of  breaking  hydrogen  bonds  in  the 
base  pair  (~3  kcal/mol). 

The  data  uncover  that  a  significant  portion  of  UDG’s 
intrinsic  binding  energy  is  used  to  drive  the  unfavorable 
process  of  DNA  bending.  Thus,  a  key  question  is  the  benefit 


4  This  is  similar  to  our  previous  conclusion  that  enthalpically 
destabilizing  the  uracil  base  pair  promotes  enhanced  binding  by 
increasing  the  probability  of  extrahelical  states  that  are  recognized  by 
the  enzyme  (14).  However,  our  previous  interpretation  did  not  include 
the  effect  of  base  pair  flexibility  on  the  energetics  of  DNA  bending. 
Indeed,  increased  flexibility  appears  to  be  more  important  than  the 
enthalpic  benefit  of  removing  hydrogen  bonds. 
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Figure  8:  Free  energy  diagrams  depicting  the  increasing  base  flipping  equilibrium  in  each  intermediate  along  the  reaction  coordinate  for 
uracil  flipping  and  the  effect  of  base  pair  destabilization  on  binding  affinity  of  UDG.  (A)  Free  energy  diagrams  depicting  the  increasing 
equilibrium  constant  for  uracil  flipping  along  the  reaction  coordinate.  The  equilibrium,  ATfiip  =  [open]/[closed],  describes  uracil  in  a  closed, 
hydrogen-bonded  state  in  the  DNA  duplex  and  an  open  state  that  is  further  along  the  reaction  coordinate  for  base  flipping.  This  flipping 
equilibrium  is  qualitatively  depicted  for  the  free  DNA,  the  two  enzyme-bound  intermediates  in  Figure  IB  (Ii,  I2),  and  the  fully  extrahelical 
product  of  the  flipping  reaction  (P).  The  known  structural  and  dynamic  features  of  each  depicted  species  are  listed  below  each  profile.  The 
equilibrium  constant  ATnjp  has  been  measured  for  the  free  DNA  (~10-5)  and  for  DNA  bound  in  Ij  (~1(T3)  using  NMR  solvent  magnetization 
transfer  methods  at  10  °C  (7),  establishing  that  UDG  increases  the  equilibrium  by  about  100-fold  in  this  early  intermediate.  Less  is  known 
about  the  stmcture  of  I2,  but  the  present  results  indicate  DNA  bending  and  largely  intact  base  pair  hydrogen  bonds.  In  the  P  state,  the  uracil 
is  held  tightly  in  the  active  site,  and  the  equilibrium  is  pushed  all  the  way  to  the  fully  open  state  (180°  rotation).  Access  to  the  closed  state 
requires  passage  back  over  the  transition  state  connecting  I2  and  P.  The  estimated  rotation  angles  for  the  closed  and  open  states  are  based 
on  the  pseudodihedral  reaction  coordinate  defined  by  Banavali  and  MacKerell  (11, 14),  which  ranges  from  10°  for  uracil  in  B  DNA  to  180° 
for  full  rotation  of  the  base  out  of  the  major  groove.  (B)  Free  energy  effect  of  base  pair  destabilization  on  binding  affinity  of  UDG  (AAGbind 
=  AAGnex  +  AAGHbond).  The  relative  free  energy  profiles  for  the  stable  UF-D  duplex  and  the  flexible  UF-M  duplex  are  depicted.  The 
increased  flexibility  of  the  UF-M  substrate  decreases  the  free  energy  of  I2  because  less  binding  energy  is  required  to  bend  the  DNA  at  this 
step.  The  magnitude  of  the  flexibility  benefit  is  estimated  as  AAGnex  ~  3  kcal/mol  (see  text).  In  addition  to  the  flexibility  benefit,  the 
absence  of  base  pair  hydrogen  bonds  in  UF-M  as  compared  to  the  three  hydrogen  bonds  in  UF-D  enhances  binding  of  the  P  state  by 
AAGHbond  ~  3  kcal/mol  (see  Figure  3B  and  text).  This  benefit  arises  because  for  UF-M  there  is  no  energetic  price  for  breaking  these  bonds 
in  the  base  pair.  For  simplicity  we  have  shown  that  AAGn'x  arises  entirely  at  the  I2  state,  but  this  effect  is  more  likely  realized  at  both  the 
I2  and  P  states. 


derived  from  this  significant  energetic  expenditure.  Since  the 
rate  of  spontaneous  base  pair  opening  is  already  fast  com¬ 
pared  to  events  on  the  base  flipping  pathway,  it  is  mecha¬ 
nistically  reasonable  to  propose  that  the  role  of  enzyme- 
induced  DNA  bending  is  to  promote  forward  commitment 
along  the  pathway  rather  than  to  lower  activation  barriers  to 
flipping  (i.e.,  stabilization  of  extrahelical  intermediates).5 


UDG  likely  promotes  forward  migration  of  the  uracil  down 
the  pathway  by  forming  nonspecific  interactions  with  the 


5  Computational  studies  on  UDG  and  other  enzymes  have  suggested 
that  DNA  bending  facilitates  base  flipping  (9-12, 51—55),  and  a  recent 
computational  study  on  Hhal  (cytosine-C5)-methyltransferase  indicates 
that  this  enzyme  acts  by  forming  a  web  of  interactions  with  several 
high-energy  intermediates  along  the  flipping  pathway  (56). 
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DNA  backbone  which  become  progressively  stronger  during 
the  formation  of  two  discrete  intermediates  (Ii  and  h,  Figure 
IB).  In  the  earliest  nonspecific  intermediate  (Ii),  NMR  imino 
exchange  experiments  indicate  that  UDG  increases  the 
opening  equilibrium  ([open]/[closed])  by  almost  100-fold  (2.8 
kcal/mol),  largely  by  decreasing  the  closing  rate  (7).  In  the 
subsequent  intermediate  (I2),  the  magnitude  of  the  equilib¬ 
rium  is  currently  unknown,  but  a  further  increase  is  suggested 
because  the  uracil  has  substantially  decreased  stacking 
interactions  with  the  adjacent  2-AP  base  arising  from  DNA 
bending.  Since  the  presence  of  uracil  is  required  to  achieve 
I2,  then  this  intermediate  must  be  considered  a  specific 
complex  in  which  the  enzyme  recognizes  some  feature  of 
the  uracil  base.  The  extent  or  nature  of  the  discrimination 
between  uracil  and  other  normal  DNA  bases  at  this  inter¬ 
mediate  step  is  currently  not  known.  In  conclusion,  our  view 
based  on  these  findings  is  that  UDG  facilitates  base  flipping 
by  utilizing  binding  energy  to  alter  the  DNA  structure, 
thereby  progressively  shifting  the  flipping  equilibrium  toward 
the  open  state  that  is  ultimately  trapped  by  specific  interac¬ 
tions  located  within  the  uracil  binding  pocket.  The  specific 
interactions  are  formed  only  very  late  in  the  process,  and 
thus,  the  earlier  steps  in  the  flipping  pathway  are  largely 
driven  by  nonspecific  interactions.  A  further  experimental 
goal  with  UDG  is  to  unambiguously  identify  the  amino  acid 
side  chains  that  promote  the  stepwise  transition  from  the 
closed  to  open  state. 

Conclusion.  We  have  uncovered  the  origins  of  high-affinity 
binding  of  UDG  to  an  extrahelical  uracil  base.  By  discretely 
altering  base  pair  hydrogen  bonding  and  site  flexibility,  we 
have  found  that  high  binding  affinity  is  determined  more  by 
the  flexibility  of  the  site  than  by  the  enthalpic  benefit  of 
removing  base  pair  hydrogen  bonds.  Thus  a  significant 
contribution  to  enzymatic  base  flipping  is  the  unfavorable 
energetic  cost  of  DNA  bending.  This  cost  is  estimated  in 
two  ways:  (i)  the  stiffening  effect  of  pyrene  (3.7  kcal/mol) 
and  (ii)  the  additional  favorable  energetic  effect  on  binding 
of  removing  the  last  hydrogen  bond  from  the  uracil  base 
pair  (~3  kcal/mol).  In  the  future,  these  simple  features  of 
high-affinity  target  site  binding  by  UDG  may  be  utilized  to 
direct  UDG  to  specific  sites  in  genomic  DNA. 
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